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I,  Introduction 

This  is  the  final  report  on  the  research  program  "Resonances 
of  Radar  Targets  and  Target  Discrimination" ,  supported  by  the 
Naval  Air  Systems  Command  under  Contract  No.  M00019--79-C-0270 ,  and 
covering  the  period  March  6,  1979  to  March  5,  19£0. 

X 

The  topic  of  the  research  program  covers  methods  for  target 
discrimination  by  radar  scattering,  making  use  of  the  natural 
resonances  of  the  targets.  The  general  plan  of  the  program  consists 
in  a  study  of  resonances  of  targets  of  simple  shape  (e.g.,  sphere) 
but  complex  composition  (e.g.,  dielectric  coating  of  conductors) 
during  the  first  year,  followed  by  simple  conducting  targets  and 
*  comparison  with  parallel  measurements  performed  at  NCWC  Dahlgren  in 
the  second  year,  with  further  extension  to  conductors  o f  more  complex 
shapes,  using  Waterman's  methods,  later  on. 

The  emphasis  of  our  study  is  two-fold,  namely:  1)  to  study 
ways  how  the  resonance  structure  of  the  radar  echo  can  be  used  for 
identifying  the  nature  of  the  target  (the  "inverse  scattering  problem") 
and  2)  to  obtain  a  physical  understanding  of  the  phenomena  that  cause 
the  resonances,  and  how  they  manifest  themselves  in  both  steady- 
state  and  pulsed  echoes.  Such  an  understanding,  in  preference  to  a 
more  mechanical  use  of  the  singularity  expansion  method,  will  provide 
us  with  the  true  capability  of  mastering  the  target  discrimination 
problem. 


f*  -  - 
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II.  Progress  Achieved 

The  first  year's  research  has  concerned  steady-state 
studies  of  the  resonances  of  dielectric  spheres  and  dielectric- 
coated  conducting  spheres.  In  the  following,  superscripts  refer 
to  the  list  of  publications  which  were  prepared  under  the  auspices 
of  this  contract,  and  which  are  listed  in  Section  III. 

1)  Dielectric  Sphere^-.  The  Theory  of  Resonance  Scattering2 
developed  by  the  Principal  Investigator  and  his  collaborators  has 
been  applied  to  this  problem  in  order  to 

(a)  obtain  resonance  frequencies  (as  a  function  of  dielectric 
constant)  in  a  purely  real  calculation,  in  contrast  to  previous 
tedious  solutions  of  the  complex  problem 

(b)  interpret  the  resonances  physically  in  terms  of  the  phase 
matching  of  circumferential  waves,  and 

(c)  obtain  the  properties  (in  particular,  the  dispersion 
curves)  of  these  circumferential  waves. 

Thi3  work  has  been  submitted  for  publication  in  IEEE  Trans.  API 
and  has  been  reported3  at  the  National  Radio  Science  Meeting, 

5-8  November  1979,  Boulder,  Colorado. 

2)  Dielectric-Coated  Perfectly  Conducting  Sphere.  This 
problem  has  been  approached  from  the  stand  point  of  inverse  scattering, 
in  order  to  show  how  the  dielectric  constant  and  the  thickness  of  the 
coating  can  be  obtained  from  the  spacing  and  the  widths  of  the 
resonances*-?,  ^  complete  resonance-theoretical  treatment  has  been 
performed  on  this  topic® separating  the  scattering  amplitudes 
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into  resonant  terms  and  a  non-resonant  background,  and  obtaining 
the  dispersion  curves  of  the  cirf umferential  waves  that  propagate 
in  the  coating  and  cause  the  resonances  , 

Our  work  on  this  problem  has  been  prepared  for  publication 
in  Journ.  Appl.  Phy3^  and  other  aspects  of  it  in  Radio  Science9 , 
and  it  will  be  reported  at  the  North  American  Radio  Science/URSI 
meeting,  2-6  June,  1980,  Quebec,  Canada5,  (abstract  attached)  as 
well  as  the  Spring  Meeting  of  the  American  Physical  Society, 

28  April-1  May,  Washington,  D.  C.  (abstract  attached) .  It  will 
further  form  part  of  an  Invited  Paper  for  a  special  Issue  on 
"'Inverse  Methods  in  Electromagnetics"  (W.  M.  Boerner,  Editor)  , 

IEEE  Trans.  A-P,  Spring  1981.  A  copy  of  this  latter  paper  is 
attached  to  the  present  report. 

3)  The  Transient  Problem  of  Conducting  Bodies. 

A  collaboration  with  an  experimental  group  at  the  Naval 
Surface  Weapons  Center,  Dahlgren,  VA  led  by  Dr.  Bruce  Hollmann 
has  been  initiated,  with  pulsed  measurements  to  be  carried  out  at 
the  Dahlgren  radar  range  on  conducting  bodies  of  simple  shapes,  to 
be  compared  to  theory  (both  the  Singularity  Expansion  Method  or 
SEM,  and  Creeping  Wave  Theory).  Meanwhile,  we  obtained  the  theory 
for  an  experiment  on  a  conducting  sphere,  and  have  established  the 
connection  between  SEM  and  Creeping-Wave  theory,  with  the  latter 
providing  physical  insight  into  the  SEM  results.  Our  analytical 
results  are  now  being  programed,  and  will  be  published  for 

acoustic*0  aa  well  as  radar  targets*!  of  spherical 
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shape  in  the  near  future.  Similar  studies  axe  now  being 
undertaken  by  A.  Nagl  and  F.  Cannata  for  conducting  infinite  or 
finite  cylinders. 

III.  Personnel 

The  following  personnel  has  been  involved  in  the  work 
described  in  tne  present  progress  report: 

H.  Uberall,  Professor,  Catholic  University,  Principal  Investigator 
A.  Nagl,  Catholic  University,  Research  Associate 
P.  J.  Moser,  Catholic  University,  Graduate  Student  (also  Naval 
Research  Laboratory,  Washington,  D.  C.) 

F.  Cannata,  Catholic  University,  temporary  Research  Associate 
J.  D.  Murphy,  Annandale,  VA,  voluntary  collaborator 

L,  Flax,  Naval  Research  Laboratory,  ’Washington,  D.  C. 

G.  C.  Gaunaurd, Naval  Surface  V7eapons  Center,  White  Oak, 

Silver  Spring,  MD 

L.  R.  Dragonette,  Naval  Research  Laboratory,  Washington,  D.  C. 

J.  Severns,  Naval  Research  Laboratory,  Washington,  D.  C. 
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III  Publications  and  Abstract-;  of  Talks  Supported  or  Partly 
Supported  by  Present  Contract, 

1.  J.  D.  Murphy,  P.  J.  Moser,  A.  Nagl,  and  If.  Uberall,  "A  Surface 
Wave  Interpretation  for  the  Resonances  of  a  Dielectric  Sphere", 

IEEE  Trans.  AP  (accepted  for  publication) 

2.  L.  Flax,  G.  C.  Gaunaurd,  and  H.  Uberall,  "Theory  of  Resonance 
Scattering",  Physical  Acoustics  Vol.  15,  Academic  Press, 

New  York,  1980  (in  press) . 

3.  J.  D.  Murphy,  P.  J.  Moser,  A.  Nagl,  and  H.  Uberall,  "A  Surface 
Wave  Interpretation  for  the  Resonances  of  a  Dielectric  Sphere", 
paper  B4-3,  International  Union  of  Radio  Science,  National 
Radio  Science  Meeting,  5-8  November  1979,  Doulder,  Colorado. 

4.  G.  Gaunaurd,  H.  Uberall,  and  P,  J,  Moser,  "Resonances  of 
Dielectrically- Coated  Conducting  Spheres  and  the  Inverse 
Scattering  Problem",  to  be  published  in  Journal  of  Appl.  Phys. 

5.  P.  J.  Moser,  J,  D.  Murphy,  A.  Nagl,  H.  Uberall,  and 

G.  C.  Gaunaurd,  "Resonances  of  a  Dielectrically  Coated 
Conducting  Sphere:  Surface  Waves  and  the  Inverse  Scattering 
Problem",  paper  submitted  to  the  North  American  Radio  Science 
Meeting  (URSI,  Commission  B) ,  Quebec,  Canada,  2-6  June,  1980. 

6.  G.  C.  Gaunaurd  and  II.  Uberall,  "Solution  of  Inverse  Scattering 
Problems  in  the  Resonance  Case",  paper  submitted  to  the 
Washington  Meeting  of  the  American  Physical  Society, 

28  April  -  1  May,  1980. 

7.  G.  Gaunaurd,  II.  Uberall,  and  L.  R.  Dragonette,  "Solution  of  the 
Inverse  Electromagnetic,  Acoustic,  and  Elastic  Scattering  Problem 
in  the  Resonance  Case",  Invited  paper  for  the  IEEE  Transactions  AP, 
Special  issue  on  "Inverse  Methods  in  Electromagentics"  (W.M.Boemer 
Editor) ,  Spring  1981. 
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8.  J.  D.  Murphy,  A.  Nagl,  and  II.  Uberall,  "Resonance  Effects 
in  Radar  Scattering  from  Conducting  Dielectric,  and  Coated 
Targets",  Navair  Electronic  Research  Program  Review,  Stanford, 
University,  Stanford,  CA,  6  March,  1979. 

9.  P.  J.  Moser,  J.  D.  Murphy,  and  H.  Uberall,  "Resonances  and 
Surface  Waves  in  Conducting  Spheres  with  Dielectric  Coating" , 
Radio  Science  (to  be  published) . 

10.  H.  Uberall,  and  J.  D.  Murphy,  "Acoustic  Surface  Wave  Pulses 
and  the  Ringing  of  Resonances",  J.  Acoust.  Soc.  Am.  (to  be 
prepared  for  publication) . 

11.  H.  Uberall  and  J.  D.  Murphy,  "The  Physical  content 

of  the  Singularity  Expansion  Method",  IEEE  Trans.  A-P  (to  be 
prepared  for  publication) . 
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Washington,  DC  20064 
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The  theory  of  resonance  scattering,  developed  by  some  of 
the  present  authors  in  the  context  of  the  scattering  of 
acoustic  waves  from  elastic  obstacles,  or  of  elastic  waves 
from  fluid-filled  cavities  or  solid  inclusions,  is  here 
applied  to  the  problem  of  radar  scattering  from  dielectrically- 
coated,  perfectly  conducting’ spherical  targets.  The 
numerous,  sharp  and  narrow  resonances  obtained  in  previous 
calculations  of  the  corresponding  radar  cross  sections  (see, 
e.g.  ,  Ruck  et  al.  Radar  Cross  Section  Handbook,  Plenum 
New  York,  1970)  are  shown  here  to  correspond  to  the  real 
resonance  frequencies  of  the  target,  which  can  be  calculated 
from  a  real  characteristic  equation  and  agree  approximately 
with  the  real  part  of  the  natural  frequencies.  A  surface- 
wave  interpretation  of  these  resonances  is  made  entirely  on 
the  basis  of  real  analysis,  thus  avoiding  the  complexities  of- 
the  Watson  transformation.  The  dispersion  curves  for  these 
surface  waves  are  obtained,  which  show  (for  the  TM  mode)  a 
discontinuous  transition  from  conductor-type  at  low  frequencies 
to  dielectric-type  at  high  frequencies.  The  inverse  problem 
is  also  solved  (see,  e.g.,  G.  C.  Gaunaurd,  H.  {Jberall,  and 
L.  R.  Dragonette,  IEEE  Trans.  Antennas  Propag.,  Special  Issue 
on  Inverse  Methods  in  Electromagnetics,  W.  M.  Bocrner,  ed., 
to  be  published)  by  showing  how  from  the  spacing  and  the 
widths  of  the  resonances,  the  dielectric  constant  and  the 
thickness  of  the  coating  may  be  obtained.  Finally,  the 
relationship  between  the  present  resonance  theory,  and  the 
singularity  expansion  method  or  SEM  (see,  C.  E.  Baum, 
Interaction  Note  88,  December  1971)  is  briefly  sketched. 
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Abstract 


Conventional  methods  of  attacking  the  inverse  scattering 
problem  are  either  based  on  the  geometrical  optics  limit,  or  on 
the  U3e  of  Fredholm  integral  equations  of  the  first  kind.  We 
demonstrate  here  that  for  the  (quite  common)  case  that  resonances 
are  present  in  the  echoes  of  electromagnetic,  acoustic  or  elastic 
waves  scattered  from  corresponding  tarots,  the  theory  of  resonance 
scattering  recently  developed  by  us  represents  a  powerful  tool  for 
the  solution  of  the  inverse  scattering  problem.  For  example,  for 
the  case  of  elastic-wave  scattering  from  fluid-filled  spherical 
cavities  in  solids,  the  spacing  of  successive  resonances  determines 
the  sound  speed  of  the  fluid-filler,  and  the  resonance  widths  its 
density.  Measurements  of  the  resonances  were  performed  by  us  for 
acoustic  scattering  from  submerged  elastic  objects.  For  the  case 
of  electromagnetic  scattering  from  dielectric  or  coated  conducting 
targets,  spacing  and  widths  of  the  resonances  are  shown  to  determine 
the  dielectric  constant  and  the  coating  thickness  of  the  target 
object  in  a  direct  fashion. 
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I.  Introduction 

Inverse  scattering  is  of  importance  in  many  fields  of  Physics, 

In  quantum  mechanics,  inverse  scattering  methods  try  to  determine 
the  potential  if  the  scattering  cross  section  is  known.  Geophysicists 
attempt  to  determine  the  distribution  of  density  and  of  seismic 
velocity  with  depth,  from  observed  reflection  responses.  In  general, 
inverse  scattering  methods  try  to  obtain  information  about  a  source 
distribution,  a  scattering  obstacle  or  inhomogeneity  of  a  medium 
from  remotely-sensed  field  data. 

Attempts  at  a  solution  of  the  inverse  scattering  problem  have 
conventionally  proceeded  via  the  use  of  the  geometrical-optics 
limit,  or  else  via  an  exact  treatment  .that  leads  to  Fredholm 
integral  equations  of  the  first  kind^-.  in  the  present  review,  we 
wish  to  demonstrate  that  if  resonances  are  present  in  the 
scattering  cross  sections  as  a  function  of  frequency,  then  the 
measurable  characteristics  of  these  resonances  furnish  us  with  a 
powerful  tool  for  the  solution  of  the  inverse  scattering  problem, 
i.e.  for  the  determination  of  the  geometry  and  the  composition  of 
the  target  from  the  measured  properties  of  the  echo  returns.  In 
fact,  the  presence  of  resonances  in  scattering  echoes  seems  to  be 
the  rule  rather  that  the  exception,  as  confirmed  for  the  cases  of 
acoustic2  and  elastic-wave-^  f  ag  well  as  electromagnetic* 
scattering. 

The  theory  of  resonance  scattering,  as  developed  by  us 
recently^  8^  provi^es  a  formalism  that  may  be  used  for  a 
straightforward  and  direct  solution  of  the  inverse  scattering  problem. 


Although  this  will  be  demonstrated  in  the  following  on  the  basis 
of  a  few  simple  examples  only  (reflection  of  acoustic  waves  from 
plane  layers;  of  elastic  waves  from  fluid-filled  spherical  cavities 
and  of  electromagnetic  waves  from  dielectrically-coated  conducting 
spheres) ,  it  is  clear  that  an  elaboration  of  these  procedures  for 
the  case  of  more  complex  targets  will,  with  proper  attention  to 
detail,  provide  us  with  a  wealth  of  information  on  both  the 
geometry  and  the  composition  of  the  target  which  hitherto  seems  to 
have  been  largely  disregarded.  Our  examples  show  that,  e.g.,  the 
spacing  of  successive  resonances  in  elastic-wave  echoes  directly 
determines  the  sound  speed  in  the  fluid  filler  of  a  spherical 
cavity,  and  that  the  resonance  widths  determine  the  density  of  the 
fluid;  or  that  for  the  case  of  electromagnetic  scattering  from 
spherical  dielectric  or  coated  conducting  targets,  the  spacing  and 
widths  of  the  resonances  determine  the  dielectric  constant  and  the 
coating  thickness  of  the  target  object  in  a  direct  fashion.  The 
resonances  thus  characterize  the  scattering  object  as  if  they  were 
its  signature,  and  their  spectrum  constitutes  a  code  identifying 
the  target  in  a  way  that  resembles  optical  spectroscopy. 

II,  Acoustic  Resonances 

The  first  example  to  be  discussed  is  the  acoustic  transmission 
and  reflection  for  a  plane  fluid10  or  elastic11  layer  of  thickness 
A  imbedded  in  an  ambient  fluid.  The  latter  will  be  characterized 
by  its  density  f  and  sound  velocity  c,  while  the  (elastic)  layer 
is  assumed  to  have  density  ,  and  compressional  and  shear  speeds 
and  t  respectively,  A  plane  sound  wave  be  incident  at  an 


angle  0  with  the  normal,  and  be  refracted  at  angles  o(  for 

compressional,  or  /3  for  shear  waves  info  the  layer.  The 
transmission  (T)  and  reflection  coefficients  (R)  may  be  written  in 
the  form 


<la) 


(lb) 


where 


(j>/C./C05  & 


_JEL  ~  si*  9)** 

?A  cos  0 


(lc) 


/h.  =  'C  ,  'YL^  -  /sc±  . 


(Id) 


For  a  lossless  layer,  one  has  jTlz  +  jRj1  -  l.  The ("symmetric* ,  s, 
and  "antisymmetric",  a)  coefficients  Cs  and  Ca  are  given  in  the 
literature**,  for  a  fluid  layer  0  , /J->0  > ,  one  finds 


C1,  -»  J , 


(2a) 


' ta.y\  £ 


where 


(2b) 


using  the  abbreviations  X  =  kd  =*(od/c,  y  =*■  sin  0  .  If  n<n.<l, 

A  t 

the  critical  angles  of  total  reflection 


0^  *  St>rX  01^  , 


0^  =  SL-w 


n, 


(2c) 


appear. 

Equations  (la,  b)  show  that  each  of  the  two  terms  in  T  as 
well  as  in  R  will  exhibit  a  resonance  behavior  when  the  real  parts 
of  the  denominators  vanish* 


or 


(3a) 


These  equations  also  describe  the  free-plate  dispersion  relations10'11 
(corresponding  to  no  fluid-loading,  r«*  0)  for  the  characteristic 
waves  on  the  plate.  The  resonances  will  be  the  more  pronounced,  the 
smaller  the  imaginary  part  X  of  the  denominator  (rouohly  speaking, 

<<  D#  i.e.  the  lighter  the  fluid  loading.  For  the  fluid 
layer,  Eqs.  (3a)  have  the  solutions 


-  7 


cot  5  =  0,  3  =  ^  =  (ws*  i  >7r 


(3b) 


or 


-t<m  3  «  0  , 


^  ^  , 


(3c) 


with  mg  and  ma  being  integers. 

The  two  independent  variables  of  the  problem  may  be  taken 
as  x  and  y.  Resonance  forms  of  T  and  R  may  be  obtained  in  either 
of  these  variables,  holding  the  other  one  fixed. 

Taking  y  fixed,  we  expand  the  functions  Cg  or  Ca  around  their 
resonance  values: 

cs,a(*)=  tx"xms,a^ cs,ax  (xm g,a)' 

where  Xms^a  designate  the  solutions  of  Cg^a(x)  “  0,  and 
cs,ax  “  ^Cs^a/9x.  For  the  reflection  coefficient,  e.g. ,  this 
leads  to  the  resonance  form 


» 


(5a) 


with  each  resonance  expression  being  valid  in  the  vicinity  of  its 
own  resonance  position  x  -  x  For  the  elastic  plate,  Eqs.  pa) 

have  to  be  solved  numerically  for  xm  ,  for  the  fluid  plate,  one 
has 


x 


5,A 


2  3  ( yi  ,l  -  stV  £  ) 

in,  ,  ** 


(5b) 


The  resonance  widths  are 


+ 


2t  /  ) 

•  S,«t  *  v  ' 


(5c) 


which  for  the  fluid  plate  become  independent  of  mg^a: 

r  =  L*p(p4  c«se)- 

These  resonances  may  be  interpreted  as  poles  in  the  lower  half  of 
the  complex  x-  plane,  located  at 


'W. 


-  i 


r 


»/<* 


X 


*  X 


(5c) 
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Similar  resonance  expressions  may  be  derived  for  R,  as  well  as 
for  T  and  R  as  functions  of  the  y-  variable. 

Equation  (5a)  i3  illustrated  in  Fig.  1  for  a  v.’ater  layer  in 
alcohol,  at  various  values  of  ®  .  Figure  2  illustrates  resonances 
in  the  angular  variable  of  the  transmission  loss  TL  =  log  j T  \2 
for  a  steel  plate  in  water.  The  figures  show  both  the  results  for 
the  individual  resonance  terms,  Eq.  (5a),  and  as  obtained  from  the 
exact  expression,  Eq.  (la) .  The  value  of  the  resonance  formalism 
lies  in  the  fact  that  it  furnishes  explicit  expressions  for  the 
position,  width,  and  strength  of  each  resonance  contribution, 

Eqs.  (5b)- (5d) ,  as  a  function  of  all  the  material  parameters. 

It  thus  appears  possible  to  carry  out  an  inverse  determination 
of  the  material  constants  of  a  plate  by  examining  the  resonances 
experimentally.  In  fact,  for  a  fluid  layer,  Eqs.  (5b)  and  (3b,c) 
show  that  the  resonance  positions  measure,  via  nd  *  c/cd,  the 
sound  speed  cd  in  the  layer,  and  Eq.  (5d)  shows  that  the  width  of 
the  resonances  measures  the  density  ratio  f  /9AjL  .  This  is  the 
simplest  example  for  the  solution  of  the  inverse  scattering  problem 
via  a  study  of  the  resonances.  For  the  elastic  plate,  the  inverse¬ 
scattering  solution  proceeds  analogously  and  is  only  slightly  more 
involved.  While  for  the  fluid  layer,  Eq.  (5b)  is  physically 
understood  by  the  kinematic  condition  that  the  normal  component  of 
the  layer  wave  number  equals  k<|  *  mTT/d,  corresponding  to  an  integer 
number  of  half-wavelengths  across  the  layer,  this  pciture  is  no  longer 
true  for  elastic  layers  where  shear  waves  enter  also.  Here,  the 
physical  understanding  is  provided  by  the  coincidence  of  the  trace 
velocity  of  the  incident  wave  with  the  speed  of  the  characteristic 
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waves  propagating  along  the  free  plato^.  However,  Eqs,  (5c) 
or  (5d)  which  conplete  the  inverse  problem,  are  only  obtained 
from  a  dynamic  calculation. 


Ill  Elastic  Wave  Resonances 

In  this  section,  we  shall  discus 3  an  application  of  the 
Resonance  Theory  of  (visco) elastic  wave-scattering  that  solves 
an  inverse  scattering  problem  of  elastodynamics ,  The  procedure, 
like  that  discussed  in  the  preceding  section,  set3  the  pace  for 
extensions  to  more  general  situations,  as  well  as  for  direct 
transfer  to  other  neighboring  scientific  areas. 

The  amplitudes  of  backscattered  waves  returned  by  obstacles 
in  viscoelastic  solids  ,  when  plotted  as  a  function  of  freqequency, 
exhibit  so  many  rapid  oscillations  and  complicated  features  that 
until  very  recently  it  wa3  not  possible  to  extract  the  physical 
information  contained  in  them  about  the  nature  of  the  obstacle. 

We  can  analyze  these  backscattered  "echoes"  using  the  Resonance 
Theory  of  viscoelastic  wave-scattering  from  cavities  in  solids,  and 
obtain  from  them,  for  a  given  shape  of  the  cavity,  the  material 
composition  of  the  filler  substance  (assumed  fluid  in  this  example, 
for  simplicity) . 

When  a  (spherical)  cavity  filler  i3  sent  into  oscillation  by 
elastic  (say,  compressional ,  or  p-type)  incident  waves,  a  3et  of 
modal  resonances  (fundamental  and  overtones)  are  excited  in  it,  that 
characterize  the  filler  a3  if  it  were  its  signature.  The  way  these 
resonances  can  be  U3ed  for  material  discrimination  purposes  resembles 
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the  way  chemical  elements  are  identified  from  their  optical  spectra. 

The  far-field  backscattering  amplitude  of  p-waves,  returned  by 
a  spherical  fluid-filled  cavity  of  radius  a  in  a  solid  under  p-wave 

n 

incidence,  is  given  by  the  "partial  wave  sum" 


i />>-  £ h>> -  r  <-*>■£*  (c-i) 

/»V=<9  'A-J  A 


(6) 


where  p-wave  velocity)  is  the 

compressional  wave  number  in  the  ambient  solid,  which  is  complex 
if  the  30 lid  is  viscoelastic  (-4#  C)  .  The  S-matrix  element  snPP 

is  given  in  the  literature7.  Figure  3  shows  the  plot  of  the 
modulus  of  the  summed  amplitude  in  the  backscattering  direction 

Q  =*  IT  (i.e.,  |^PP(7T  )/a|  )  for  a  cavity  filled  with  ethy  1-alcohol 
in  an  aluminum  matrix.  This  is  the  "echo"  containing  the  rapid 
oscillations  and  complex  features  mentioned  above.  Fifteen  partial 
waves  were  added  to  produce  that  plot,  which  for  the  ran qe 
0  4  a  4  10  of  Fig.  3,  yields  enough  accuracy.  This  is  the 
amplitude  as  a  sensor  would  record  it. 

An  analysis  of  the  first  partihl-wave  contributions  making  up 
the  sum  is  shown  in  Fig.  4.  Each  mode  can  be  split  into  a 
background  (center  graph)  and  resonances  (bottom  graph) .  An  index 
X  labels  each  resonance  of  each  mode  n.  Overlaying  the  resonances 
which  we  ju3t  isolated,  on  top  of  the  summed  amplitude  of  Fig.  3, 
identifies  each  extremum  in  that  plot  with  the  sets  of  modal  reso- 
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nances  contained  in  the  partial  waves.  The  resonance  responsible 
for  each  wiggle  of  Fig.  3  is  labeled  by  indices  (n,  h  ,  and  the 
graph  shews  over  thirty  wiggles  uniquely  identified  with  nodal 
resonances  in  this  fashion.  The  location  of  all  the  resonance 
spikes  in  fact,  identifies  the  ratio  c^/c^  of  the  f iller-to-r.atrix 
wavespeeds.  The  filler's  sound-speed  cr  can  be  found  from  the 
spacing  A  between  any  two  high-order  consecutive  overtones 
shown  in  Fig.  4.  For  spherical  cavities  in  solids,  we  have  shown^-2-14 
that  for  £»  i,  one  obtains 


A  a  i r  (C//^)  . 


(7) 


For  alcohol  in  aluminum,  that  relation  gives  A=  0.59,  just  as 
observed  in  Fig.  4.  Thus,  knowing  cd  for  the  matrix,  and  the 
spacing  A  between  consecutive  high-order  modal  resonances, 
determines  the  sound  speed  Cf  of  the  filler,  and  vice  versa. 

A3  in  the  acoustic  case,  the  resonance  theory  of  scattering 
lead 3  to  a  resonance  form  of  the  scattering  amplitude7 


(8) 


s  • 


13  - 


where  Sn^°^^  3  exp  (2-t^-n.)  represents  the  S-  matrix  element  of  a 
"soft"  (i.e.,  evacuated)  cavity.  The  first  term  in  brackets  o* 

Eq.  (8  )  leads  to  the  resonances  (bottom  portion  of  Fig.  4)  ,  and 
the  second  term  to  the  non-resonant  background  (center  portion  of 
Fig.  4) ;  their  coherent  addition  produces  the  interference  effects 
in  each  partial  wave  (top  portion  of  Fig.  4) .  The  resonance 
frequencies  *n£  refer  to  a  normalized  frequency  scale  x  =  a . 

The  density  ratio  of  f iller-to-matrix  materials  (i.e.,  p/fv  > 
is  found  from  the  resonance  widths.  Explicit  forms  for  the  widths 

f~~p  of  the  resonances  may  be  given? »  12-14 
'n<-  w  *  and  were  snown  to 

depend  on  S?  /  ^  f  ,  via  the  mechanical  impedance  Fn(x)  of  the 

obstacle  which  is  a  known  quantity.  If  we  evaluate  the  latter  at 

a  point  one  half-width  below  any  resonance  peak  (i.e., 

x  “  xn£  “  /2J  and  tlien  expand  for  x  >>  1,  the  result  is 


- 


-f- 


7T‘ 


&x 


■+ 


(9) 


(  n ,  X?  *  integers  )  f 

which  holds  for  integer  n  and  ,  so  that  consecutive 

resonances  do  not  overlap.  In  the  high-frequency  limit,  this  gives 
the  filler-to-matrix  density  ratio  in  terms  of  the  uniform  asymptotic 

spacing  b  ,  and  the  width  of  any  (high-index)  overtone  X  ,  of 

any  mode  n,  which  are  either  all  known  or  previously  determined 
quantities . 
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3  5 

For  ethyl  alcohol  (  =  0.79  g/cm  ,  cf  =  1.213  x  10  c n/sec) 

3  5 

in  aluminum  (^  =  2.7  g/cm  ,  =  6.420  x  10  cm/sec)  we  find 

A  —  0.59 ,  and  for  the  first  (i.e,  n  =  1)  mode,  the  Graphical  width 

of  its  tenth  overtone  {.£.  =  10}  ,  which  occurs  at  x^  =  £A  =  5.93, 

is  approximately  0,06,  as  found  from  a  plot  sinilar  to  the  bottom 

portion  of  Fig.  4.  To  obtain  the  width  as  defined  in 

12  r~ 

Eq.  (8),  we  must  still  divide  by  f3,  and  substitution  into 
Eq.  (9)  yields  ^  / <^>  =  0.296,  which  contains  a  negligible  error 
when  compared  to  the  true  ratio  0.79/2.70.  The  (fluid)  filler  is 
completely  identified  once  its  sound  speed  and  density  are 
determined  by  this  asymptotic  procedure.  We  conclude  that  the 
sound-speed  ratio  between  filler  and  matrix  is  extracted  from  the 
spacing  between  any  two  consecutive  overtones  of  any  mode ,  and  the 
matrix- to-fi Her  density  ratio  is  extracted  from  the  width  o^  any 
high-order  resonance  by  means  of  Eq.  (9).  Thus,  the  high-frequency 
region  X  #1,  £>>1  contains  all  the  information  about  the  material 
composition  of  the  filler.  The  filler  has  been  identified  from  its 
"echo",  and  an  inverse  scattering  problem  has  been  solved  in  a  novel 
fashion . 

Echo  spectra  from  fluid  filled  cavities,  such  as  that  shown  in 

Fig.  3,  have  been  measured  by  Pao  and  Sachse  using  ultrasonic 

techniques.  In  acoustics  ,  measurements  of  echoes  from  submerged 

elastic  objects  have  been  performed  using  short-pulse  experimental 
2 

techniques  .  Figure  5  show3  the  (total)  far-field  scattering 
amplitude  modulus  for  an  aluminum  cylinder  in  water,  plotted  V3,  k a 
where  k  **  CO/c  is  tie  wave  number  of  sound  in  water,  and  a  the 
cylinder  radius.  The  measured  points  agree  well  with  an  exact 
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calculation  (solid  curve) ,  demonstrating  the  accuracy  with  which 
the  interference  minima  caused  by  the  resonances  may  be 
experimentally  determined  in  the  acoustic  case,  and  hence  the 
accuracy  with  which  the  inverse  problem  may  be  solved  here. 

IV  Electromagnetic  Uave  Resonances 

A  procedure  analogous  to  the  one  used  above  for  elastic-wave 

scattering  will  now  be  employed  in  order  to  solve  the  inverse 

scattering  problem  for  electromagnetic  waves  from  a  dielectrically 

coated,  perfectly  conducting  sphere.  For  a  plane  wave  with 

propagation  constant  k  =  CJ/c  incident  along  the  z-axis  on  the 

coated  sphere,  the  far-field  amplitude  of  the  scattered  electric 

4 

field  i3  given  by 

[  ee  -  2*^(6)  ■* ''*>*]  , 

using  spherical  coordinates,  and  designating 


1G 


The  "hie  coefficients"  are  given  by 


x^(x)  -  <  ( yj*  (x))' 


tx)  -  <  Xv  ( x  ^  ( x)) 1 
X  ^  J  }(xj  -  i'X  (x^1'  oo) 


(12b) 


where  x  -  ka,  a  being  the  (outer)  radius  of  the  spherical  target. 
For  a  perfectly  conducting  sphere,  the  modal  impedances  and 
admittances  Z^  and  Yn ,  respectively,  become  Zn  -*  O  and  ■*  <*>  . 

^"or  a  coated  conducting  sphere  (radius  b  of  the  conducting  core, 
thickness  3  m  a  — b)  which  will  be  considered  belov/,  expressions  for 
Zn  and  ^n.  are  9iven  the  literature4. 

The  eigenvibrations  of  the  target  correspond  to  the  zeroes  of 

the  denominators  of  an  and  b^,  The  resulting  characteristic 
equations 


(xC<*>)  /( x Aj1'  (X)),  (TE  nodes) 


C13a) 


(TM  modes) 


(13b) 
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lead  to  complex  eigenf requencies.  Hov.-cvcr ,  since  the  resonances 
of  the  problem  are  sharp  and  narrov?  (as  seen  from  the  following 
figures),  indicating  small  imaginary  parts  of  the  denominators, 
it  is  sufficient  to  search  for  the  zeroes  of  the  real  parts  of  the 
denominators.  The  corresponding  equations 


1  j. tZ*  =  "Re\  (TE  nodes) 

iYn,—  Re  ^  }  (T:i  nodes) 

■re 

are  solved  by  the  real  eigenfrequencies  x  =  x  . 

TM 

and  of  the  TE  and  TM  modes,  respectively. 

Our  resonance  theory  may  again  be  applied  to  the  present 
15 

scattering  problem  ,  by  writing  the  Mie  coefficients  in  the  form 


(14a) 


(14b) 


(15a) 


(15b) 


^  -  - 
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where 


u> 

(x-Hl'’  (x‘))/ 


(10) 


and  where  we  designated 


.  v-TE 


-£JV>  = 

(M 


(xiL,v’oo)'  =  ^,2-T™ 


(x  2J1,1  («)) 


(17a) 


(17b) 


these  latter  expressions  representing  the  scattering  amplitudes 
from  a  conducting  sphere.  While  thus  the  second  terms  in  brackets 
of  Eqs,  (15)  will  provide  a  non-resonant  background  to  the  scattering 
amplitudes,  the  first  terms  provide  resonances  as  in  the  previous 
sections.  In  fact,  expanding  their  denominators  about  the  resonance 
frequencies  and  ,  one  may  bring  the  Mie  coefficients 

into  the  form 


< 

i 

\ 

i 


I 


-  1*  - 


'VI 


*  •*T£‘ 

±  O 

2.  ^ 


if 


M. 


Tfc' 


x-/£^rTE 

^  *  *«e  2-  nfi 


(13a) 
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TM 


-T* 


.  .Th  . 
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exhibiting  the  resonances  directly.  From  the  mentioned  expansions, 

r-  re  r—  Tn 

explicit  expressions  for  the  widths  /  ,  I 

/viE  ♦iv 
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may  also  be  obtained  . 

The  electromagnetic  resonances  are  illustrated  in  Fig.  6  for 
the  example  of  a  perfectly  conducting  sphere,  coated  with  a 
dielectric  of  relative  coating  thickness  =*  J/a  =  0,1,  and 
dielectric  constant  £ ^  =  6.  The  backscattcring  cross  section4, 

O'-  TT  I  £  (-1)"  , 

*  '  =  4 

is  plotted  vs.  ka  up  to  25#.  numerous  resonances  are  visible  here, 
and  have  been  identified  a3  TE  or  TU  resonances  by  solving  Eqs, 
(14a)  or  (14b).  These  solutions,  denoted  by  (nE^)  or  (n  M  £ )  , 
respectively,  are  marked  by  arrows  in  Fig.  6,  thus  indicating  how 
the  resonances  may  be  read  off  a  measured  back3cattering  cross 
section  curve  obtained  a3  function  of  frequency,  in  order  to  carry 
out  an  inverse  scattering  analysis.  In  fact.  Fig.  6  may  be  viewed 


(18c) 


1 
i 


l 


***** 


as  the  "characteristic  spectrum"  identifying  the  scattering 
object. 

V7e  have  also  evaluated  numerically  some  of  the  partial  -rave 
coefficients  XX^  and  i'.\ ich  exhibit  tlie  (ni:£)  and  (nil/) 

resonances  for  a  given  partial  wave  number  n.  Figure  7  presents 
the  TE  moduli  [(2m.+  l)/n(n  +  1)3  I  •f*/*  l  for  n  =  9  (top)  and  10 
(bottom)  plotted  vs.  ka.  The  fundamental  (£=  1)  partial  v/ave 
resonance  is  clearly  visible  on  the  left  of  each  figure,  while 
the  remaining  broad  structures  are  dominated  by  the  conducting- 
sphere  background  [second  term  in  square  brackets  of  Eqs.  (18)]  , 
in  which  the  overtone  (£^2)  resonances  (first  term  in  square 
brackets)  are  imbedded.  The  senaration  of  these  resonances  from 
the  background  will  be  described  in  a  forthcoming  paper16. 

Using  the  expressions  of  Eqs.  (18) ,  and  the  known  forms  of 
Zn  andY^,  the  inverse  problem  may  now  be  solved  as  in  the 
previous  cases.  For  the  higher-order  resonances  located  in  the 
high-frequency  spectral  region  (x  >>  1) ,  one  may  show  using 
asymptotic  forms  of  the  Bessel  functions  that 


±/<  co-t  kj_  $ 


(19  a) 


l  -  VF*  ki  <5 


(19b) 
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(where  k^  *  <*Vc^  k,  c^  beina  the  speed  oe  liaht  in  the 

dielectric) ,  and  that  in  the  sane  limit,  the  right-hand  sides  of 
the  characteristic  equations  Eqs.  (14)  vanish.  Accordingly ,  their 
solutions  are 

=  (2i*i)V*  (20a) 

<20b> 

witli^  =  integer.  Ue  therefore  have  for  the  spacing  of  two 
adjacent  (TE  or  Til)  resonances  in  the  high-frequency  spectral 
region  (on  the  k-  scale) : 


(kJ) 


TE,TM 


-  (u)7'™  -wvf* 


>1 


X 


(21a) 


which  represents  one  equation  for  the  two  unknowns  cT  and  6 


of  the  inverse  problem.  Calling  ^ »  a-b  =  a(l-y),  and 

re.TM 

the  spacing  on  the  x-  scale,  Eq.  (21a)  may 


re.rn 
be  written  as 
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(21b) 


A  second  equation  for  the  unknowns  nay  be  found  from  a 
consideration  of  the  resonance  widths.  These  are  qiven  by  the 
resonance  theory  in  the  form  (e.g.  for  TE) : 


|  r 

/Vi 


T£* 


J>n  f  <«(*>) 
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When  evaluated  at  one  half-width  below  the  resonance, 

Tf  r-  TE-  v  T£ 

*  =  S'*-  —  Ag  >  anc^  using  asynptoti^  forms 

obtains  the  equation 


i.e.  at 
,  one 


=  -Lan  [  (t-f) ],  (22b) 

which  together  with  Eq.  (21b)  suffices  to  determine  the  two  unknowns 
^1  and  X  tbe  Prob^em  from  the  measured  resonance  spacings 
Ag  and  resonance  widths 
resonance  spectrum. 


Q 


in  the  high-frequency  region  of 


the 


V.  Conclusion 


The  foregoing  discussion  has  demonstrated  on  the  basis  of  soma 
selected  examples  from  acoustic,  elastic-wave  and  electromagnetic 
scattering,  that  the  presenca  of  resonances  in  scattering  cross 
sections,  and  their  analysis  by  the  resonance  theory  of  scattering 
as  recently  developed  by  us ,  may  serve  to  solve  the  inverse 
scattering  problem  in  a  rather  direct  fashion.  For  the  simple 
examples  considered  here,  all  the  unknown  material  parameters  of  the 
target  could  be  obtained  from  a  measurement  of  the  asymptotic 
spacing  and  of  the  width  of  the  resonances.  A  known  shape  (flat 
layers,  or  spherical  objects)  has  been  assumed  for  the  targets  of 
tile  examples  considered.  It  is  clear,  however,  that  the  resonance 
spectrum,  and  the  widths  as  well  as  the  intensities  of  the 
resonance  lines,  will  depend  equally  well  on  the  geometrical  shape 
as  on  the  material  composition  of  the  target  (in  fact,  for 
impenetrable  targets  such  a3  acoustically  rigid  objects  or  perfect 
conductors,  the  spectrum  will  depend  on  shape  only) .  A  complete 
experimental  determination  of  the  properties  of  all  the  resonances 
appearing  in  the  echo  returns  from  scattering  objects  should 
therefore  identify  the  object  as  far  as  its  geometrical  shape  and 
material  composition  is  concerned,  in  the  same  fashion  in  which 
optical  spectra  identify  chemical  substances  that  emit  these 
spectra.  An  entire  scientific  discipline  of  acoustic,  elastic-wave 
or  electromagnetic  spectral  analysis  could  therefore  be  developed 
for  more  or  less  complex  scattering  objects,  with  the  above-discussed 
basic  examples  to  be  considered  as  a  starting  point. 
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rigure  Captions 


Fig.  1.  Plots  of  1 T \ ^  vs  fd  at  various  angles  of  incidence  0 
for  the  exact  expr jssion  (solid  curves),  and  for  the 
individual  resonances  in  the  resonance  approximation 
(broken  curves) ;  case  of  a  water  layer  in  alcohol. 

Primed  numbers  corresponding  to  symmetric,  unpriraed 
to  antisymmetric  resonances.  Here  =  41.8102°; 
f  **  frequency. 

Fig.  2.  Transmission  loss  in  decibels  versus  $  at  73.1 

kHz-in.;  (top  portion) s  superposition  of  resonance 
forms;  (bottom  oortion) j  calculated  from  exact  theory; 
for  a  steel  plate  in  water. 

Fig.  3,  Modulus  of  summed  p-*p  backscattered  amplitude  plotted 
vs.  k^a  (adding  15  terms  in  the  mode  series) for  an 
alcohol-filled  cavity  in  aluminum.  Some  of  the  rapid 
oscillations  in  this  plot  are  accounted  for  (as  indicated) 
by  the  resonances  shown  in  the  following  figure. 

Fig.  4.  The  zeroth  (i.e.  n  **  0)  modal  contribution  (top)  of  the 

/a  |  amplitude  of  an  alcohol-filled  spherical  cavity 
in  an  aluminum  matrix,  and  its  separation  into  a  smooth 
background  (center)  and  isolated  resonances  (bottom) .  The 
resonances  in  this  n«*0  mode  are  labeled  and  marked.  The 
asymptotic  spacing  i3  seen  to  be  A  0.59. 

Fig.  5.  Comparison  of  theory  (solid  curve)  and  experimental 


i 


observation  (points)  for  the  modulus  of  the  far-field 
acoustic  back scattering  amplitude  of  an  aluminum  cylinder 
in  water. 


Fig.  6.  Backscattering  cross  section  of  a  perfectly  conducting 
sphere  coated  with  a  lossless  dielectric  of  dielectric 


constant  E^  =  6,  and  relative  coating  thickness 

=  0.1.  as  a  function  of  k  .  Electric  (nE^)  and 

a 

magnetic  (nll^ )  resonances  are  identified  by  arrows. 
Fig.  7.  Partial-wave  TE-  moduli  £(2n+l) /n(n+lf|  |  |  plotted 

vs.  ka,  for  n  -9  (top  and  10  (bottom) . 
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